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Abstract 

We analyze how the radiation pressure interaction between a mechanical el- 
ement and an intensely driven optical cavity mode can be exploited for gen- 
erating squeezed light. We study in particular how the performance of the 
optomechanical device can be improved when a homodyne-based feedback 
loop is added to control the motion of the mechanical element of the system. 
We show that, when driving the cavity at resonance, an appropriate propor- 
tional feedback control is able to improve the generation of ponderomotive 
squeezing, which should be detectable with state-of-the-art apparatuses. 
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1. Introduction 

Optomechanical sensors are extensively used for detecting small displace- 
ments, masses and forces [H, Q, S]- The recent progress in nanofabrication 
techniques now enables these devices to enter a regime where they can man- 
ifest quantum effects. As a consequence, a new and highly active research 
field has emerged with the aim of designing and implementing devices in 
which a strong and tunable optomechanical interaction allows to manipulate 
at the quantum level the state of micromechanical resonators and optical 
modes. Various quantum effects have been predicted [3j, and current efforts 
focus on optically cooling the resonator to its quantum ground state 0, 0|, 
and to the more challenging goal of generating entanglement between optical 
and mechanical modes p, |7|, |8j. 

However, the first prediction of manifestly quantum effects in cavity op- 
tomechanical system was made more than fifteen years ago 0, [To| and con- 
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cerned ponderomotive squeezing, i.e., the possibility to generate quadrature- 
squeezed light at the cavity output due to the radiation pressure interaction 
of the cavity mode with a vibrating resonator. In fact, the mechanical ele- 
ment of the cavity is shifted proportionally to the intracavity intensity, and 
consequently the optical path inside the cavity depends upon such inten- 
sity. Therefore the optomechanical system is equivalent to a cavity filled 
with a nonlinear Kerr medium, which is able to produce squeezed light at 
the output [11]. This analogy between the ponderomotive interaction and 
Kerr nonlinearity was first pointed out in 12L 1 1 31] and used to demonstrate 



radiation-pressure- induced optical bistability [12j. In an ideal Kerr medium 
the response to the optical field is instantaneous and therefore significative 
squeezing is produced with a bandwidth limited only by the cavity and with 
no excess noise. In an optomechanical device, the dynamics of the mechanical 
element becomes relevant and significantly modifies the squeezing spectrum; 
the thermal noise acting on the mechanical element represents a limitation, 
but significant squeezing can still be obtained if radiation pressure effects 
predominate over thermal noise, which requires large enough cavity finesse 
and mechanical quality factor. 

The first analysis of Refs. jsl, [lo[, based on a Fabry- Perot cavity and a 
description of the mechanical motion as a single harmonic oscillator, was 



then later extended to the case of many vibrational modes in [14|. The 
problem has been recently reconsidered in a Michelson interferometer setup 
in [l5j], and a preliminary experimental study of the possible signatures of 
ponderomotive squeezing in a Fabry- Perot cavity with a movable end-mirror 



has been recently carried out in [16 



The experimental demonstration of ponderomotive squeezing would be 
relevant both at fundamental and at a practical level. In fact, on one hand 
it would demonstrate the possibility to manipulate in the quantum regime 
microresonators with microgram masses, and on the other hand it would 
represent a new source of squeezed light, alternative to o ptic al parametric 
amplifiers, which could be used, for example, in spatial [17| and spectro- 



scopic measurements [18( . In particular squeezing is useful for improving the 



sensitivity of gravitational wave (GW) interferometers. As first proposed by 



Caves [19|, the signal to shot-noise ratio can be improved without increas- 
ing the driving power by injecting squeezed vacuum states of light into the 
signal output port. Later it was realized that injection of squeezed light 
in the interferometer can also be used to reduce the overall quantum noise 
including radiation pressure noise, thereby beating the standard-quantum- 
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limit 20|, |2l|, |22[. In order to be useful for GW interferometers, a source 
of squeezed light must work at sideband frequencies within the audio band 
(10 Hz-10 KHz), and should have an adjustable phase of squeezing (i.e., the 
phase with sub-shot noise variance). Standard parametric amplifier sources 
squeeze light at higher sideband frequencies and at a fixed phase, but recent 
works have shown an impressive advance in the adaptation of these sources 
to GW interferometers. In fact, squeezed states have been demonstrated 

and were tested on a suspended GW prototype 
Furthermore, it has been shown that detuned Fabry- 



at audio frequencies [23 



interferometer 24 



Perot cavities can convert a squeezed vacuum with frequency-independent 
squeeze quadrature into one with frequency- dependent squeeze quadrature 



25] . These filters have been shown to be broadly applicable to existing inter- 
ferometer configurations 25|, |26] and therefore provide the required tunability 
of the phase of squeezing. 

In this paper we study if and how the generation of ponderomotive squeez- 
ing can be improved by adding a suitable feedback control to the cavity op- 
tomechanical system. In fact, feedback controls have been already proposed 
and successfully implemented in these systems. A notable example is cold 
damping feedback [2J7|,l28|,l29|,|30| in which a feedback loop adds a viscous force 



able to cool the resonator, even to its quantum ground state [3 if . and which 



has been experimentally applied in various laboratories [32], l33l. l34j. l35l l36l. 137 
Moreover, feedback control has been also proposed for suppressing radiation 
pressure noise in GW interferometers through a "quantum locking" scheme 
[38| . and it is therefore interesting to verify if feedback can also help the gener- 
ation of ponderomotive squeezing. Our analysis will include all fundamental 
quantum noise sources and thermal noise, and will neglect "technical" noise 
sources, such as laser and electronic noise, even though they can be not neg- 
ligible in the audio band. We shall see that a proportional feedback control 
is able to improve squeezing in an appreciable way. 

In Sec. II we provide a quantum Langevin description of optomechanical 
systems, while in Sec. Ill we show how to add feedback controls to them. 
In Sec. IV we determine the spectrum of squeezing of the output light and 
we determine its general properties. In Sec. V we specialize to the resonant 
case and we present some numerical results, while Sec. VI is for concluding 
remarks. 
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2. Quantum Langevin description of the optomechanical system 

In a typical cavity optomechanical system one has a cavity mode inter- 
acting with a mechanical oscillator with a term which is proportional to the 
light intensity. The standard situation corresponds to a Fabry-Perot cavity 
with one heavy, fixed mirror through which a laser of frequency ujq drives 
a cavity mode, and another light end-mirror of mass m (typically in the 
micro or nanogram range), which is free to oscillate. The optomechanical 
interaction is provided by the radiation pressure of the cavity mode on the 
oscillating light mirror. However the same description applies also to related 
systems, with different geometries and in which the coupling may be pro- 
vided by different mechanisms. Notable examples are silica toroidal optical 
microcavities which are coupled to radial vibrational modes of the support- 
ing structure or partially transparent SiN membranes oscillating within 



a standard high-finesse cavity [39|, |40| . In these two systems the coupling is 



still provided by radiation pressure, while in other cases the coupling has a 



different origin, as for example in [4l|, |42|, |43[ where it stems directly from 
the dipole gradient force. 

For clarity we shall refer from now on to the prototypal situation of an 
optical Fabry-Perot cavity of length L formed by a rigid massive mirror at 
one end and a vibrating micromechanical mirror at the opposite end, even 
though the description can be easily adapted to the other configurations. 
The laser significantly drives only a single cavity mode with frequency u c , 
from which it is detuned by A = oo c — u . The motion of the micro- mirror 
can be described by the set of its vibrational normal modes, each with its 
own resonance frequency Uj and damping rate 7-,-. The Hamiltonian of the 
system can be written as 

H = hw c a^a + ^(P? + <§) + H int + ihE^e^ - ae luJot ), (1) 

3 

where the cavity field annihilation operator a satisfies the commutation re- 
lation [a, at] = 1, and the mechanical modes are described by dimensionless 
position and momentum operators satisfying [qk,Pj] = iSkj- Denoting by k 
the cavity decay rate, the parameter E is related to the input power P in by 
E = ^/2Vi n K/Hujo. Phase and amplitude laser noise could be included by 
assuming that E is a complex stochastic process with fluctuating modulus 
and phase. The single cavity mode description is valid in the adiabatic limit 
when all the relevant mechanical frequencies Uj are much smaller than the 
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cavity free spectral range c/2L, which is typically satisfied for short cavities. 
In this limit the scattering of photons by the mirror motion from the driven 



mode to other cavity modes is negligible 44]. The interaction between the 
cavity mode and the vibrational modes is described by H int and it is due to 
the radiation pressure acting on the surface S of the vibrating mirror. One 



has 14 



d 2 

'S 

where P(r) is the radiation pressure field and 



H mt = - / d 2 rP(r)-u(r), (2) 



u 



= ( 3 ) 



is the displacement field of the mirror surface at point r. This field can be 
written as a sum over the corresponding (dimensionless) displacement field 
of each normal mode, itj(r), which is characterized by an effective mass rrij = 
p J d 3 r |-Uj(f)| 2 (p the mirror mass density). We consider a one- dimensional 
situation, i.e., we assume that the driving laser and the cavity are perfectly 
aligned. In this case, light is sensitive only to mirror surface deformations 
along the cavity axis, u x (r), so that Eq. (j2J) becomes 



H int = - / d 2 rP x (r)u x (r). (4) 
Js 

In general, the radiation pressure due to an optical power V impinging on a 
mirror with reflection coefficient 1Z can be written as 

P.if) = -^Kv 2 opt (r), (5) 

with v op t(f) denoting the spatial structure of the incident optical field on the 
mirror surface. Within the cavity, one can rewrite TPjc = h(u c /L)a'a and 
also assume 1Z ~ 1. One ends up with 

H^t = -hy]G 3 a}aqj, (6) 

3 

where the optomechanical couplings are given by 



Gi = ^J— , (7) 
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and 

Cj = / (frvl^iu^^r) (8) 
Js 

is the overlap at the mirror surface between the cavity mode and the j-th 
mechanical mode. Due to the chosen normalization of v% t (r) and Uj(f), 
the overlaps satisfy the condition —1 < Cj < 1. Eqs. ©-0 show that 
the radiation pressure directly couples the cavity mode only with the mir- 
ror collective displacement operator q e ff = J2j^o1j- When the detection 
bandwidth involves only a single, isolated, vibrational normal mode of the 
microresonator, the collective coordinate q e ff is well approximated by the 
selected normal mode, and the single harmonic oscillator description usually 
adopted is justified. In the most general case, one has to include in the dy- 
namical description of the system all the vibrational normal modes which 
contribute to the detected signal. 

The unavoidable action of damping and noise onto the dynamics associ- 
ated with the Hamiltonian of Eq. ([TJ) is described by adopting the formalism 



of quantum Langevin equations |45|, |46( which, in the frame rotating at the 
laser frequency uq, are given by 

Qj = UjPj, (9) 
Pj = -Ujlj ~ l&j + GWa + £j, (10) 
a = -(K + iA )a + i^2G J a qj + E + V2^a in . (11) 

3 

The input noise a in (t) describes the optical vacuum field entering the cavity, 
and it is delta correlated in the time domain {a in (t)a m ^(t')) = 5(t — t') 
[III, while the mechanical Brownian stochastic forces with zero mean value 
are uncorrelated from each other and have the following, generally non- 



Markovian, correlation functions |46 



< oth ( \ 1 

1 2k B T 1 



(12) 



with ks the Boltzmann constant and T is the temperature of the reservoir 
of the micromechanical mirror. 

We are interested in ponderomotive squeezing and therefore one requires 
a strong radiation pressure interaction, which is achieved when the intra- 
cavity field is very intense, i.e., for high-finesse cavities and enough driving 
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power. In this limit (and if the system is stable) the system is character- 
ized by a semiclassical steady state with the cavity mode in a coherent state 
with amplitude a s (\ct s \ ^> 1), and a new equilibrium position for the vibra- 
tional modes, displaced by g|. The parameters a s and q 3 s are the solutions of 
the nonlinear algebraic equations obtained by factorizing Eqs. (j9"|)- (fTT|) and 
setting the time derivatives to zero. They are given by 



r<3 



a: 



|2 



qi = -^. t (13) 



A = 0, (14) 

E , . 

a s = — r-rr, (15) 
K + zA 

where the effective detuning A is obtained from A by subtracting the fre- 
quency shift caused by the steady state radiation pressure 

A = A -|a s | 2 V^E. (16) 

j J 

These steady state equations are responsible for the optical bistability ob- 



served in [12J and analyzed in [13j. Then, we linearize Eqs. (19|)- (|TT|) around 
the steady state values by writing operators as sums of averages plus fluctu- 
ations: a = a s + 5a, qj = q J s + 5qj and pj = p 3 s + 5pj. The nonlinear terms 
5a) 5a and 5a5qj can be ignored when the fluctuations are much smaller than 
the mean value, and this is certainly satisfied when \a s \ ^> 1. One therefore 
arrives at a system of linearized quantum Langevin equations 

5 qj = UjSpj, (17) 

5pj = -UjSqj - 'jjSpj + GjSX a + £j, (18) 

5X a = -K5X a + A5Y a + V2iiX?, (19) 

5Y a = —K,5Y a — A5X a + ^ Gj5qj + V2kY™. (20) 

We have chosen the phase reference of the cavity field so that a s is real and 
positive, we have defined the cavity field quadratures SX a = (5a + 5a))/\/2 
and 8Y a = (5a — 5a))/i\/2 and the corresponding Hermitian input noise 
quadratures X* n = (a in + a in *)/y/2 and F a in = (a in - a in ^)/iy/2. We have 
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also defined the effective optomechanical couplings 



Gj = G{a s %f2 = ^p-J Vi ;1 \ o\ ' (21) 

3 L y m^-wo + A 2 ) v 7 

3. Adding the feedback loop 

The fluctuation dynamics described by Eqs. (IT71)-( l20l) is then modified by 
adding a feedback force acting on the vibrational modes of the resonator (see 
Fig. 1), that can be applied in various ways, either exploiting the radiation 
pressure force of an additional laser beam (as in Ref. |32|), or through elec- 



tromechanical actuators as in |34j, |35|, |36|, |37J • The feedback loop is obtained 



by extracting a fraction of the cavity output which is then processed in order 
to drive an appropriate actuator acting on the resonator. The simplest and 
most efficient way to extract the feedback loop mode is by means of a beam 
splitter with amplitude transmission t and reflection r, with t 2 + r 2 = 1 (we 
can always choose mode phases so that t and r are real), so that 

d = ta out - rb in , (22) 
c = ra out + tb in . (23) 

d is the annihilation operator describing the overall output of the device, 
and c is the annihilation operator of the light employed for the feedback 
loop. This latter optical mode is detected by a balanced homodyne detector 
measuring with quantum efficiency rj the quadrature with phase 9. Detection 
with non-unit efficiency rj can always be described in terms of an effective 
beam splitter with amplitude transmission ^frj in front of a perfect detector 



45] , so that the detected mode is given by 



S = y/r]C + a/1 - TjVin, (24) 

where Vi n is the annihilation operator describing the optical vacuum noise 
unavoidably entering the detector. In practice the bosonic mode s describes 
the photocurrent at the output of the homodyne detector which is then elec- 
tronically filtered and amplified in order to actuate the vibrational modes 
(see Fig. 1). 

The feedback loop generally acts differently on the different mechanical 
modes and it is described by an additional force term on Eq. ( TTSj) . given by 



8 




Figure 1: Schematic description of a cavity optomechanical system subject to a homodyne- 
mediated feedback control of the mechanical resonator. 

the time convolution 

{«*)} /6 = - f dt'g^t-t'f-jf,, (25) 

where gj(t) is the causal feedback transfer function on the j-th vibrational 
normal mode and 9 s (t) = (se~ ld + s^e ld ) / a/2 is the detected field quadrature. 
The homodyne phase 9 and the transfer functions gj are the feedback loop 
parameters which must be optimized in order to achieve the best possible 
squeezing of the output mode d. From E qs. (|2"2"| )-( ]2"4"|) and using the input- 
output relation a out (t) = y/2na(t) — a,i n (t) [45], one has 

s(t) = y f 2KTjra(t) - y/rjra in (t) + y/rjtb in (t) + a/1 - rjv in (t). (26) 
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Inserting Eq. into Eq. §Z5§ and then adding the resulting feedback force 
into the quantum Langevin equations of Eqs. ([17]) - (1201) . we arrive at the 
following Fourier-transformed equations 

— iu5qj(u) = U)j5pj(u)), (27) 
— iu5pj((jj) = —ujj5qj((jj) — jj5pj(oj) + [Gj — ry/rjgj((jj) cos#] 5X a {oj) 

-r^jg 3 (u) Bm96Y a (uj) + + ^ln fb (uj), (28) 

V 2k 

-iw<Df a (a;) = -«£X (w) + A<^ a (w) + V2^ n (w), (29) 

-iw<jy«(a;) = -K<5r tt (w) - A6X a (u) + ^ G^w) + v / 2^F a in (w),(30) 

i 

where 

= M n (") - vWM - v / i 3 ^CH (si) 

is the vacuum input noise injected by the feedback loop, and we have defined 
the Fourier transform of the field quadratures of a given mode with annihi- 
lation operator f(u), 6f{uj) = [f(u)e- i6 + p(u)e i6 ]/V2. Notice that, due to 
the fact that we have written the quantum Langevin equations in the frame 
rotating at the laser frequency u , the frequency u is referred to this latter 
frequency, i.e., u = corresponds to the laser frequency. 

The observables of interest are the Fourier transform of the field quadra- 
tures of the output mode d given by Eq. ( 122]) . The explicit expression of 
these quadratures as a function of the cavity, mechanical, and feedback loop 
parameters can be obtained by replacing into Eq. ([22]) the input-output re- 
lation and the explicit solution of Eqs. ( 127]) - (|30]) for the cavity field ampli- 
tude quadrature 5X a {u) and the phase quadrature 5Y a (cu). After long, but 
straightforward calculations one gets the Fourier transform of these quadra- 
tures as a linear combination of all the noise terms acting on the system 

X d (u) = ^(w)!;^) + a 2 ( W )C(w) + <T3(w)ir(w) + <T4(w)n m (w) 

+* B (u)e*(u) + ae(u)J2Gixf\uM"), (32) 

j 

Y d (u) = m(u)X^(u) + fi 2 (u)Y^ n (uj) + n 3 {uj)Xi n {uj) + ^(u)Y* n (u) 

+/x 5 (u;)CM + toWY^GjxfMtji"), (33) 

j 

where the coefficients are given by 

ai(u)) = tD(w) -1 [(k + iu) [k - icj + ry/r] sin 9 \ g (u)] - A [A - A G (w)]p4) 
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a 2 {uj) = tADiuy 1 [2k + r-y/rj sin 9Xg(oj)} 

+r [A 2 -AA G H]], 
= -Dicu)- 1 [t 2 Ay^jsm9\g(cu 



(35) 

[Ay 7 ?/ cos 8 X g (co) + r (k — ico) [k — iu + r^/rj sin 8X g (co)] 

(36) 

0-4(0;) = -D[u)~ l \r&y/7iwn9\g{u)\ , (37) 
(T 5 H = -tADiu^^/T^Xgiu), (38) 
cr 6 (u) = tAD(u)~ 1 V2K, (39) 

for the output amplitude quadrature and 



/i 2 (w) 
/i 4 (o;) 



-tD(o;) _1 [2kA - 2kA g + (k + iu)r y/fj cos 0\ g {u)] 



-D(uj)-H 2 [(« - 1^)^/77 cos 0A ff (u;)] , 

— Z^(w) -1 [(k — ico)y/rj sin 8X g {u) 

+r [(« - iw) 2 + A 2 - AA G (w) + Arv^cos0A ff (w)]] 



7x5(0;) = —tD(u) 1 (k — 20;) a/ 1 — TjX g {uj) 
fj,e(u) = tD(u)~ L {k — iuj)\/2k 



(40) 



tD(uj)- 1 [k 2 + o; 2 - A [A - A G (w) + rv^cos0A a (w)]] , (41) 



(42) 

(43) 
(44) 

(45) 



for the phase quadrature. We have used the definition 



D{uj) — (k — ico) [k — ico + ry / r/sin^A s (o;)]+A [A — A g (oj) + r^fr] cos 8 X g (uo)} . 

(46) 

and we have introduced the two frequency-dependent quantities 



X G (co) = ^G 2 X fM, 



(47) 
(48) 



depending upon the unperturbed susceptibility of the j-th vibrational mode 



CO: 



co 2 — co 2 — iuj^j 



(49) 



The function X 9 {uj) describes the effect of feedback, while A G (o;) plays the 
role of an effective mechanical susceptibility of the resonator (see also Q). 
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4. The squeezing spectrum 

In order to characterize the present optomechanical device as a source 
of squeezed light we have to determine the quadrature noise spectrum of 
the output mode d which, for a generic quadrature <f>d(w) is defined, due to 
stationarity, by the relation 

(M")M<J)) = St(co)S(co + u'). (50) 

The output light is squeezed at phase <p when S%(u) < 1/2, i.e., when the 
noise spectrum is below the shot-noise limit (equal to 1/2 with our defini- 
tions). The quadrature noise spectrum S^(u) can be written in terms of the 
amplitude and phase quadrature noise spectra, Sf(uj) and Sj(u) respec- 
tively, and their correlation spectrum S' ( f y (w)0, as 

s» = SM±jIM + sJM^M cos ^ + ^ H sin H . (51) 

However, rather than looking at the noise spectrum at a fixed phase of the 
field, one can perform an optimization and consider, for every frequency u>, 
the field phase 4> op t(u) possessing the minimum noise spectrum S%(oj): in this 
way one defines the optimal squeezing spectrum which can be easily obtained 
by minimizing Eq. fl5Tj) . 



Sf(u) = mm S*(u) = \ \s?(lj) + Sj{u) - yJ[S*(u,) - Sj(u)] 2 + 4 [S™(u,)] 

2SX(co)Sj(u)-2[Sr 



(52) 



Sf(u>) + Sj(u) + ^[S^u)-Sj(u)] 2 + 4[S^(u)] 2 

The frequency- dependent optimal phase is correspondingly given by 

1 



^opt(w) = - arctan 



2 ™ (53) 



We have to remember that these noise spectra satisfy the Heisenberg uncer- 
tainty theorem, originating from the fact that the two quadratures X d and 



l Sf y (lj)S(lj+lj') = [{X d (u)Y d (w')) + (Y d (u')X d (u)) + {X d (u')Y d (w)) + (Y d (w)X d (u>'))] /4 
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Yd are non commuting observables. This fact constraints their variances and 
the corresponding noise spectra according to [48[ 

S^u)Sj(u)-[S^ Y (u)} 2 >\. (54) 

The squeezing properties of the device are fully determined by the three noise 
spectra S*(u), Sj(cu), and S^i^u), which can be obtained from Eqs. (|32|) - 
(|33|) and the knowledge of the spectrum of all the noises affecting the system. 
Since the various noise terms are all uncorrelated, one has the following 
general expressions 

S*(cj) = Mu,)| 2 SJ» + \a 2 {u)\ 2 Sl m (u) + 2Re KH^u;)} S™(cj) 

+ I^HI 2 Sg» + |a 4 (u;)| 2 + 2Re {a s (o;)oj(w)} S$£(u) 

+ \a 5 (u)\ 2 S 6 v>in (u) + K( W )| 2 coth (J^) Im{A G (w)} , (55) 
5j(w) = |/ii(w)| 2 Sj» + | M2 (u;)| 2 Sjjw) + 2Re S$£(u>) 



+ \^(u)\ 2 Sl m (u) + l/ieHrcoth (^^J Im{A G (a;)}, (56) 
Y ( W ) = Re W^KM} S*» + Re {a 2 (u,K(u,)} 

+Re{ ( r 1 (a;)^(a;) + ffJ M/»!M}S2H + Re^H^H) ^n(^) (57) 
+Re {a 4 (w)/4(u;)} S 6 y in (u;) + Re {a 3 (u)nl{u) + <7 4 (w)/i5(w)} S£ Y {u) 

+Re{a 5 (u;)n*{uj)}S e v ^(uj) + Re KM/i*(u;)} coth (-^) Im{A G (u;)}. 

Here we have assumed generic spectra for the input noise entering the cavity 
a in and for the one entering the beam splitter b in . We have also employed the 
thermal Brownian noise correlation function of Eq. ( TT2|) and the definitions 
of Eqs. ( 148|) - (|49|) . The fact that the thermal Brownian noise spectrum is 
proportional to Im{A G (a;)} is a consequence of the fluctuation-dissipation 
theorem (47| . 

These equations give the most general expression for the output spectrum 
of squeezing for a cavity optomechanical system subject to feedback. The 
best ponderomotive squeezing is achieved by minimizing S° d pt (u) of Eq. (|52}. 
However, such a minimization is subject to two constraints. The first con- 
straint comes from stability conditions. In fact, we are considering the quan- 
tum fluctuations of the optomechanical system around its stationary state, 
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which is possible only if the system parameters are within the stability re- 
gion. When the system is stable there is no exponential divergence in its 
time evolution and this is possible only if all the poles of the solutions in 
the frequency domain of Eqs. ( I32|) - (l33j) have negative imaginary part. From 
the explicit solutions of the coefficients <Ji(u) and one easily sees that 

stability is guaranteed when all the zeros of D{uj) of Eq. (|46p have negative 
imaginary part. Eq. (|46|) shows that these stability conditions are generally 
quite involved and strongly depend upon the explicit form of A s (o;), i.e., of 
the feedback transfer functions gj{oo). The second constraint comes from the 
fact that all the feedback transfer functions gj(t) must be causal, as it is 
evident in Eq. (1251) . 

5. Ponder omotive squeezing at resonance 

Let us now apply the general description of the preceding sections to the 
case when the driving laser is locked at resonance with the cavity mode, i.e., 
when A = 0. First of all we consider the standard situation of input optical 
noises a in and b in and v in in the vacuum state, implying 



Sa,in( U ) = S lin(u) = S^u) = S£m(w) = S£ in (w) = y ( 58 ) 



The resonant case is particularly convenient not only because the description 
considerably simplifies, but especially because the limitations imposed by the 
stability conditions are much less stringent. In fact in this case the relevant 
zeros are the solution of 



The sign of the imaginary part of the solutions depends upon the explicit 
form of Xg(u}), but it can be verified that stability is much easier verified 
with respect to the off-resonant case. In particular the system is always 
stable without feedback, r\ g (u) = 0, or when sin# = 0. 

The most relevant simplification of the resonant case is that the amplitude 
fluctuations 5X of the cavity field are decoupled from the other variables 
(i.e., both from phase fluctuations and from the resonator motion); as a 
consequence one has simply 



1 




(59) 




(60) 




(61) 
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as it can be verified from Eq. ( 15 5 p and the fact that when A = 0, 02(0;) = 
cr 4 (o;) = (f^ius) = a 6 (cu) = 0, <Ti(u) = t(tz + ilo)/{k — icj), and a 3 (co) = 
— r. When S d (u) = 1/2 it is convenient to rewrite in a new form both 
the Heisenberg inequality of Eq. ( )54|) and S°, pt (ui). In fact, the Heisenberg 
inequality becomes Sj(u) > 2[S d XY {oj)} 2 + 1/2, which suggests the following 
parametrization 

5j(w) = 2[5'J y H] 2 + i + 5S(w), (62) 

so that the Heisenberg condition becomes simply S d (u) > with S d (cj) a 
"residual" spectrum measuring the distance from the minimum uncertainty 
condition. In the resonant case the optimal squeezing spectrum can be rewrit- 
ten as 

s o P t {uj) = S* d {u) + 1/2 

1 + S5(«) + 2[S*r( u )]2 + y/[S* d (u) + 2[S™(u)]*] 2 + 4[S™(u)] 2 ' 

(63) 

showing two important aspects of ponderomotive squeezing at resonance. 
First of all one has that the optimal squeezing spectrum is squeezed as soon 
as S d Y {uj) 7^ 0. However this result is not relevant by itself in practice 
because, as mentioned in the Introduction, the present treatment neglects 
technical noise contributions, such as electronic noise and phase and ampli- 
tude noise of the driving laser. This technical noise adds to the expected 
noise spectrum and, if S° d pt (uj) is not too much below the shot-noise limit, it 
prevents the observation of squeezing. Therefore, one has to make S d pt (oo) 
as small as possible, and Eq. (1631) suggests how this can be achieved. In 
fact, the strongest squeezing is obtained when the two limits S d (u) <C 1 and 
[S,f y (u;)] ^> 1 are simultaneously satisfied: 



* 777^77^2 for W « K°»] 2 » L (64) 
8[#d Ml 



This strong squeezing is obtained at the optimal phase 

, * 1 r 1 

4>o P t{u) ~ -- arctan , 
1 l^d \ UJ )_ 



(65) 



which is very close, but strictly different from = 0, where the corresponding 
quadrature Xa is just at the shot-noise limit (see Eq. (IBTj) ). This means that 
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at fixed frequency, squeezing is achieved only within a narrow interval for the 
homodyne phase around </> op t(co), of width 



uj) ~ 2 \(f) opt {u)) \ ~ arctan 



This extreme phase dependence is a general and well-known property of 
quantum squeezing, which is ultimately due to the Heisenberg principle: the 
width of the interval of quadrature phases with noise below the shot-noise 
limit is inversely proportional to the amount of achievable squeezing. This 
implies that in order to detect squeezing one has to tune and stabilize the 
phase of the homodyne detection apparatus with extreme accuracy. 

By using the general equations of the former section, one gets the ex- 
plicit expressions for S* y (uj) and S^(u) at resonance as a function of the 
parameters of the optomechanical system: 

= nt 2 Re\ Xg ^ . 1 , (66) 

I [K + too) \ k — too + sin uX g {uj)\ I 

cr ( ^_ 2r \ Q XY ( ^2 | 2Kt 2 Im{A G (a;)}coth(^;/2AC B r) 
- ~W\- b d + — ; ; — ; . ,2 — 

1 \K ~ luJ + r yV sm "X g {Ul) I 



+2K 2 t 2 Im 



Ag(^) 



(k + itu) [k — itu + ^/fjr sm9\ g (u) 



t 2 


X g (u)\ 2 (K 2 + u 2 ) - 4:Kry/rj cos 6Re {\* G (uo)X g ( 


oj){k + ioj)} 


+ 2 


{k 2 + UJ 2 ) 


k — ioj + ry/rj sin 8 X g (u) 


2 



(67) 



These two equations show how to choose the system parameters in order to 
optimize ponderomotive squeezing. The fact that S^(co) must be as small as 
possible while keeping S* y (cj) very large implies first of all that one has to 
take t — > 1, i.e., the output beam splitter must possess high transmissivity, 
so that the first contribution of Eq. fl67|) to S r d (u) is negligible. This means 
that only a tiny fraction of the cavity output light is used in the feedback 
loop, but this can always be compensated by adjusting the feedback gain, 
i.e., the modulus of X g (uj). The second term of Eq. f lBT|) is the thermal noise 
contribution to the spectrum, which is smaller at low temperatures. However, 
apart from lowering temperatures, the most efficient way to suppress the 
thermal noise contribution is to choose a mechanical resonator with normal 



16 



modes possessing large quality factors Qj = Uj/jj- In fact 

Im{A G (u;)} = J2 , 2 Gj2 %T 2 2 , (68) 

so that 

G 2 

Im{A G (u;)} ~ ^$^^77"' w < (69) 

Im{A G (o;)} = u = Uj . (70) 

j J 

This means that choosing large values of Qj has the effect of concentrating 
the thermal noise contribution only within the frequency bands correspond- 
ing to the very narrow (width ~ ujj/Qj) mechanical resonance peaks, while 
the value of Im{Ac(a;)} is negligible in a wide low-frequency bandwidth 
below the mechanical resonances. Also the third contribution of Eq. f l67|) 
becomes very small when Im{Ac(a;)} is small and therefore a high-quality 
mechanical resonator is able to suppress also this contribution to the out- 
put homodyne spectrum. Therefore we expect that the optimal squeezing 
spectrum is well below the shot-noise limit in a wide low-frequency band be- 
low the mechanical resonances, while being at shot-noise level at mechanical 
resonance frequencies. 

The fourth contribution to S^(cj) is related to the optical input vacuum 
noise injected in the system by the feedback loop, and it describes the main 
effect of the feedback control on the squeezing spectrum. This fourth term 
may become negative, showing that feedback may enforce squeezing, i.e., one 
may have feedback-assisted ponderomotive squeezing. Considering frequen- 
cies u exploiting that Ag(w) is essentially real over a wide frequency 
range, and that r C 1, one has that the feedback-induced contribution to 
the spectrum is minimum (and negative) when 

\ g (u;)~2r^jcose\ G (uj), (71) 

which can be satisfied for different choices of the feedback transfer functions 
Qj{oj) in general. If we look for a solution valid at all frequencies u, the 
simplest choice is to take 

g j {co)~2r^jcos9G j Vj, (72) 
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that is a frequency-independent proportional control, which is causal and 
very different from the derivative (viscous) control employed in cold damping 



schemes for cooling the vibrational modes |32|, |34|, |35|, |36|, |37| ■ Eq. ( 172]) defines 
the optimal feedback control for generating ponderomotive squeezing, which 
is however difficult to implement because it requires gj oc Gj, Vj. In fact, this 
means a feedback transducer which couples with each vibrational mode of the 
resonator in the same way as the cavity mode does. This situation could be, 
at least approximately, realized by applying the feedback control on the me- 
chanical element through the radiation pressure of an additional laser beam 
with the same spatial transverse profile of the cavity field. However, the gen- 
eral condition of Eq. (ITTj) is easier to satisfy, without requiring the condition 
of Eq. (172 p . if one wants to improve ponderomotive squeezing only within 
a not too broad frequency interval. For example, at low frequencies, below 
all mechanical resonances u <C Uj, the unperturbed susceptibilities xf\°°) 
of Eq. f l49|) and therefore Ag(w) are frequency-independent, so that Eq. f JTTj) 
is satisfied by taking again a proportional control and simply adjusting the 
overall gain of the feedback loop. 

The noise subtraction caused by the optimal feedback of Eq. (172"]) is 

AS r d ~ -4r 2 r? cos 2 #4 ~ cos2 6 [ S d Y ] * , (73) 

which however can never be too negative; in fact, under optimal conditions, 
even though Sf Y — > oo, one has r 2 [5^' y ] — > 0, which stems from the 
fact that S^(u) must approach zero (see the first term in Eq. fj67j) ). This 
feedback-induced noise subtraction is due to destructive interference between 
the intracavity field fluctuations and the fluctuations injected by the feedback 
loop and transferred into the cavity by the optomechanical coupling. 

Feedback also modifies the cavity response as described by the expression 
of D(u) of Eq. (H6|) . and this modification of the cavity and mechanical res- 
onator response is just at the basis of the feedback cooling process proposed in 
23, 28, 29[ 30, 31] and experimentally demonstrated in 32|, 33, 34, 35, 36, 37 



On the contrary, we have seen that feedback modification of cavity response 
is not useful for the generation of squeezing and its only relevant effect is 
of making the system more easily unstable (see also Eq. fl60|) ). This means 
that the feedback optimal for squeezing corresponds to 9 = 0, because in this 
case there is no feedback-induced modification of the cavity, the system is 
always stable, and at the same time the feedback-induced noise subtraction 
of Eq. f !73|) is maximum. Choosing 8 = means measuring the amplitude 
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quadrature of the mode reflected at the output beam splitter, which at res- 
onance, is not coupled to the mechanical element, and therefore in this case 
the improvement of squeezing due to feedback and described by Eq. ( 175]) 
is not determined by the feedback modification of the cavity dynamics but 
only by the destructive interference between the intracavity fluctuations and 
those injected by the feedback loop. 

The performance of such optimized feedback control (t = 0.99, 6 = 0, 
r] = 1) is shown in Fig. 2, where the optimal squeezing spectrum (u) is 
shown (full line) and compared to the case without feedback (dashed line). 
The plot refers to a cavity with bandwidth k ~ 1 MHz, corresponding to 
a finesse F ~ 8000, length L = 6 cm, driven by a laser at 1064 nm and 
with input power P in = 30 mW. We have considered a mechanical resonator 
with a number of vibrational modes with resonance frequencies between 150 
and 600 KHz, all with the same effective mass rrij = 100 ng, and with the 
same quality factor Qj = 10 4 , placed at a temperature T = 4 K. We see 
that the optimal squeezing spectrum is well below the shot noise limit out of 
the mechanical resonance peaks, both with and without feedback, due to the 
fact that with the chosen parameter values, the radiation pressure interaction 
predominates over thermal noise. Feedback control provides a visible, even 
if not macroscopic, noise reduction over a wide frequency range. In Fig. 2b 
the same curves in the low-frequency band below the resonance peaks are 
shown, and one sees that the noise spectrum is roughly three times smaller 
due to feedback. 

In Fig. 3a we show the optimal field phase (fi op t(u) corresponding to the 
optimal squeezing spectrum shown in Fig. 2. The curves with and with- 
out feedback are indistinguishable and in both cases the optimal phase is 
practically frequency-independent in the low-frequency band up to 100 KHz. 
Fig. 3b shows the noise spectrum of a given quadrature with phase 4> (at the 
fixed frequency uj = 10 KHz), versus 0, for the same parameters of Fig. 2. 
At the small, nonzero value corresponding to 4> opt (oj) of Fig. 3a, one gets up 
to 20 dB of ponderomotive squeezing, but only within an extremely narrow 
phase interval around opi (w), as it always occurs for squeezing, since the 
width of the useful phase interval is inversely proportional to the maximum 
achievable squeezing. 
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Figure 2: Optimal spectrum of squeezing in the presence (full line) and in the absence of 
feedback (dashed line) . Feedback improves squeezing in a visible way. We have considered 
a cavity with bandwidth k = 1 MHz, length L — 6 cm, driven by a laser at 1064 nm 
and with input power P„ = 30 mW. The mechanical resonator has 25 vibrational modes 
with resonance frequencies between 150 and 600 KHz, all with the same effective mass 
Ttij = 100 ng, and with the same quality factor Qj = 10 , placed at a temperature T — 4 
K. We have also taken optimal conditions for the feedback control: t — 0.99, rj = 1, 9 = 0, 
gj ~ 2ry/rj cos 8Gj = 0.28Gj. Fig. 2b is a zoom of the low frequency region below 120 
KHz. 



6. Conclusions 

We have presented a general description of ponderomotive squeezing and 
discussed if and how an added feedback control may improve squeezing at the 
cavity output. Appreciable ponderomotive squeezing is achieved in general 
when the noise correlations caused by radiation pressure predominates over 
thermal noise, i.e., when 2[S* Y (u)] 2 ^> S^(u) (see Sec. V). Using the above 
results this condition can be written as 

2[Sj°»] 2 ^ [Re{A G (a;)}] 2 ^ V m uj 7 2 Q >> 

S r d {yS) k coth(fia;/2fc B T)Im{AG(a;)} mc 2 uj' 2 n n ' 

where in the last ratio we have assumed u ~ 0, considered for simplicity a 
single vibrational mode with mass m, quality factor Q, frequency u m , and 
we have denoted with n = kBT/hu m its mean thermal vibrational number. 
T = ttc/2kL is the cavity finesse. 

We have seen that by adding a suitable feedback control one can get a 
moderate improvement of ponderomotive squeezing thanks to a destructive 
interference between the intracavity light fluctuations and those injected by 
the feedback loop. The optimal feedback control is obtained by employing 
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Figure 3: a) Optimal field phase opt (u;) corresponding to the optimal squeezing spectrum 
shown in Fig. 2. The curves in the case with and without feedback are indistinguishable, 
b) Noise spectrum (in dB) of the quadrature at phase <f> at the fixed frequency uj = 10 KHz, 
versus the phase </>, with (full line) and without feedback (dashed line). Parameters 
are the same as in Fig. 2. 

an highly transmitting output beam splitter (r — > 0), by measuring the 
amplitude quadrature (9 = 0), and by adopting a proportional control such 
that gj(co) oc Gj. 
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